# Full text of "Parametric down-conversion from a wave-equations approach: geometry and absolute brightness"

## See other formats

o o (N in ^ ■ Oh. a- en > en m vn en o oo O X Parametric down-conversion from a wave-equations approach: geometry and absolute brightness. Morgan W. Mitchell^ ' ICFO-Institut de Ciencies Fotoniques, Mediterranean Technology Park, 08860 Castelldefels (Barcelona), Spain (Dated: 26 September 2008) Using the approach of coupled wave equations, we consider spontaneous parametric down- conversion (SPDC) in the narrow-band regime and its relationship to classical nonlinear processes such as sum-frequency generation. We find simple expressions in terms of mode overlap integrals for the absolute pair production rate into single spatial modes, and simple relationships between the efficiencies of the classical and quantum processes. The results, obtained with Green function techniques, are not specific to any geometry or nonlinear crystal. The theory is applied to both degenerate and non-degenerate SPDC. We also find a time-domain expression for the correlation function between filtered signal and idler fields. I. INTRODUCTION Spontaneous parametric down-conversion (SPDC) has become a workhorse technique for generation of pho- ton pairs and related states in quantum optics. Im- provements in both nonhnear materials [l| and down- conversion geometries have led to a steady growth in the brightness of these sources [3, y, 0, [a, M 0, @, Q ■ Apphcations of the bright sources include fundamental tests of quantum mechanics, quantum communications, quantum information processing, and quantum metrol- ogy [10! IHIj Ej [i3| • Although down-conversion sources typically have bandwidths of order 10^^ Hz, for the brightest sources even the output in a few-MHz window can be useful for experiments. This permits a new ap- plication, the interaction of down-conversion pairs with atoms, ions, or molecules. Indeed, sources for this pur- pose have been demonstrated [14| . Many modern ap- plications use single-spatial-mode collection, either for improved spatial coherence, to take advantage of fiber- based technologies, or to separate the source and target for experimental convenience. Remarkably, despite the importance of bright, single- spatial-mode sources, general methods for calculating the absolute brightness of such a source are not found in the literature. By absolute brightness, we mean the number of pairs per second that are collected, for specified beam shapes, pump power, filters and crystal characteristics. A number of calculations study the dependence of bright- ness on parameters such as beam widths or collection angles, but these typically give only relative brightness: the final results contain an unknown multiplicative con- stant [la, [la]- While useful for optimizing a given source, they are less helpful when designing new sources. A re- cent paper computes the absolute brightness for a specific geometry: gaussian beams in the thin-crystal limit [17| . In this paper, we calculate the absolute brightness for narrow-band, paraxial sources. The results are quite gen- eral, for example they apply equally well to crystals with spatial or temporal walk-off, for non-gaussian beams, etc. The Green-function approach we use is well suited to describing the temporal features of the down-conversion pairs, and we are able to predict the time correlations in a particularly simple way. To our knowledge, this method of deriving the time-correlations is also novel. Perhaps of greatest practical importance, we derive very simple relationships between the efficiency of classical parametric processes and their corresponding quantum parametric processes. For example, in any given geome- try the efficiency of sum-frequency generation and spon- taneous parametric down-conversion are proportional. This allows the use of existing classical calculations and/or experiments with classical nonlinear optics to pre- dict the brightness of quantum sources. The paper is organized as follows: In section |lTl we de- scribe briefly the variety of theoretical treatments that have been applied to parametric down-conversion, and our reasons for making a new calculation. In Section mil we describe the formalism we use, based on an ab- stract paraxial wave equation and Green function solu- tions. In Section [IV] we calculate the absolute brightness and efficiencies for non-degenerate and degenerate para- metric down-conversion and the corresponding classical processes. In section IVll we summarize the results. II. BACKGROUND The characteristics of parametric down-conversion light have been calculated in a number of different ways. Kleinman '18| used a Hamiltonian of the form H'^-l /d'xE-xiE (1) and the Fermi "golden rule" to derive emission rates as a function of frequency and angle. Zel'dovich and Klyshko [l9l | proposed to use a mode expansion and calculate pair rates treating the quantum process as a classical para- metric amplifier seeded by vacuum noise . Detailed treat- ment along these lines is given in [20|, l2l| . The problem of collection into defined spatial modes was not consid- ered, indeed the works emphasize that the total rate of emission is independent of pump focusing. After the observation of SPDC temporal correlations by Burnham and Weinberg [22], MoUow [23] described de- tectable field correlation functions (coincidence distribu- tions) in terms of source-current correlations and Green functions of the wave equation. This Heisenberg-picture calculation derived absolute brightness for multi-mode collection, e.g., for detectors of defined area at defined positions. It did not give brightness for single- mode collection, nor a connection to classical nonlinear pro- cesses. Hong and Mandcl [24| used a mode-expansion to compute correlation functions based on the Heisenberg- picture evolution and an interaction Hamiltonian of the form Hr = ^ d^x^^^E^E.E^. (2) As with MoUow's calculation, they find singles and pair detection rates, but only for multi-mode detection [30| Ghosh, et al. [2^ used the same Hamiltonian in a Schrodinger-picture description, truncating the time evolution at first order to derive a "two-photon wave- function." This last method has become the most popular description of SPDC, including work on efficient collec- tion into single spatial modes[l^, [la|. Many works along these lines are cited in reference [1^. Recently, Ling, et al. [l3] calculated the absolute emission rate based on a similar interaction Hamiltonian and a gaussian-beam mode expansion. In this way, they are able to calcu- late the absolute pair rate for non-degenerate SPDC in a uniform, thin crystal into gaussian collection modes. As described in Sect ion IVCl our calculation agrees with that of Ling et al. while also treating other crystal geometries, general beam shapes, and degenerate SPDC. Notable differences among the calculations include Heisenberg vs. Schrodinger picture and calculating in di- rect space vs. inverse space via a mode expansion. While they are of course equivalent, Heisenberg picture calcu- lations are easier to compare to classical optics, while Schrodinger picture calculations are more similar to the state representations in quantum information. As our goal is in part to connect classical and quantum efhcien- cies, we use the Heisenberg picture. Also, we note that the Schrodinger-picture "two-photon wave-function" has a particular pathology: the first-order treatment of time evolution means the Schrodinger picture state is not normalized and never contains more than two down- conversion photons. While this is not a problem for calcu- lation of relative brightness or pair distributions [29. |27| . it does prevent calculation of absolute brightness. The choice of inverse vs. real space calculation is also one of convenience: for large angles in birefringent media or detection in momentum space, plane waves are the "nat- ural" basis for the calculation. However, most bright sources use paraxial geometries and collection into de- fined spatial modes, e.g., the gaussian modes of optical fibers. In these situations, the advantages of a mode expansion disappear, while the local nature of the x'^' interaction makes real-space more "natural." Thus we opt for a real-space calculation. Our treatment of SPDC is based on coupled wave equa- tions, a standar d ap proach for multi-wave mixing in non-linear optics [2^. The calculations are done in the Heisenberg picture, so that the evolution of the quan- tum fields is exactly parallel to that of the classical fields described by nonlinear optics. This allows the re-use of well-known classical calculations such as those by Boyd and Kleinman [29|. As in the approach of MoUow, we use Green functions to describe the propagation, and find results that are not specific to any particular crys- tal or beam geometry. Unlike MoUow's calculation, we work with a paraxial wave equation (PWE). This al- lows us to simply relate the classical and quantum pro- cesses through momentum-reversal, which takes the form of complex conjugation in the PWE. We focus on narrow-band parametric down-conversion, for which the results are particularly simple. By narrow- band, we mean that the bandwidths of the pump and of the collected light are much less than the bandwidth of the SPDC process, as set by the phase-matching condi- tions. This includes recent experiments with very narrow filters [ij, but also a common configuration in SPDC, in which the down-conversion bandwidth is ^ 10 nm while the filter bandwidths are < 1 nm. III. FORMALISM A. description of propagation We are interested in the envelopes £± for forward- and backward-directed of parts of the quantum field E^^\t,yi) = (£+ exp[-f zA:z] -I- 5_ exp[— ifcz]) exp[— iwi] where k is the average wave-number and lo is the car- rier frequency. These propagate according to a paraxial wave equation V±E± = S- ■±; (3) where T>± is a differential operator and S± is a source term (later due to a x non-linearity). The formal (retarded) solution to equation ([3]) is £± {x) = £o± [x) + / d^x'g± {x- x')S± {x') (4) where x is the four- vector (t,x), £o±{x) is a solution to the source-free {S = 0) equation, and Q± are the time- forward Green functions, defined by V±g±{x;x') ^ 5^{x-x') Q±{x;x') =0 t<t'. (5) For illustration, we consider the paraxial wave equation (PWE), for which ,2 s± = UJ 9 ' ± ,(NL) (6) (7) Here V^^ is the transverse Laplacian, k = n(a;)a;/c is the wave- number, Vg = duj/dk^ is the group velocity, and -p(NL) jg |-]-^g envelope for the nonlinear polarization. We note that P± is invariant under translations of x, and that time reversal i — > — i is equivalent to direction- reversal and complex conjugation, i.e., 'D± —>■ V* . The results we obtain will be valid for any equation obeying these symmetries. In particular, the results will also ap- ply to propagation with dispersion and/or spatial walk- off, which can be included by adding other time and/or spatial derivatives to 7). From the symmetries of 'D±, it follows that the Green functions depend only on the difference x — x' , and that Q+{t,x;t' ,x') = ^l(i,x', f',x). Also, the time-backward (or "advanced") Green functions Ti±, defined by V±n±{x;x') = S'^{x-x') n±{x;x') =0 t>t' (8) oheyn±{x;x') ^gi{x',x). B. boundary and initial value problems If the value of the field is known on a plane z = Zgrc; the field downstream of that plane is £±{x) = f3z f d^x'g±{x-x')£±{x')5{z' - z,,,) (9) where Pz = ±2iA;. Similarly, if the field is known at an initial time t — to, the field later is £±ix) =Pt I d'^x'g±{x]x')£±{x')6{t' -to) (10) where Pt = 'iik/vg. Similar relationships hold for the advanced Green functions. If the field is known in some plane z — zq downstream, then £±ix) = PI I d'^x"H±{x]x')£±{x')5{z' - Zq) = P* j d^x' £±{x')5{z' ~ za)gUx';x) (11) while if the field is known at some time tf in the future, £±{x) = P*t I d'^x'n±{x]x')£±{x')5{t' -tf) = p; f d^x'£±{x')S{t' -tf)g*Jx';x) (12) C. quantization The field envelopes are operators which obey the equal- time commutation relation [£ix,t),£\x',t)]^A^^S^x'~x) (13) where A-y = y^hu!/2nngeo is a photon units scaling fac- tor and Ug = c/vg is the group index. For narrow- band fields, A^^ {^^^) describes a photon number den- sity, and VgA-^ {£^£) and VgsVgiA~^A-^ (£l£l£i£sj de- scribe single and pair fluxes. We find the unequal-time commutation relation from equation ([TI (14) [£{x),£\x')],^,.^PtAl^g{x-x') so that {Q\£{x)£\x') |0) = PtA^^g{x]x') for t > t' . For the PWE, A-'^Vg = IncEolhuj and PtAi^ = ihjj'^/c^£Q. To calculate singles rates, we will need to evaluate ex- pressions of the form {££^)- For this, a useful expression is derived in the Appendix: Equation (|A2p {£ix)£Hx')) 2hnuj^ c^eo d'^x"S{z" - Zq) xg*{x";x)g{x";x'). (15) Here zq is any plane down-stream of x and x' . D. single spatial modes A single spatial mode M± (x) is a time-independent solu- tion to the source-free wave equation X'±M±(x) = 0. Af|.(x) is the corresponding momentum- reversed solu- tion X';pM|.(x) = 0. We assume the normalization J d^x\M±{x)\'^6{z) — 1. For single-mode collection, it will be convenient to define the projection of a field £{x) onto the mode M as (16) £M{t) = / d-'xM*{y:)S{z - zo)£ix) (here and below, the +/— propagation direction is the same for £,M). Here zo is some plane of interest, and £m (t) describes the magnitude of the field component in this plane. Similarly, if the envelope is constant, the field distribution is £ix)^£Mit)M{^). (17) The optical power is (MKS units) PM{t) — 2nceo J d^x\£{t,x)\'^d{z - zq) ~ 2nceo\£M{t)\'^- Given an upstream source S{x), the M component of the generated field is £M{t) = I d^xd^x'M*{x)S{z~ zo) xg{x;x')Six'). (18) If the source is time- independent, then Equation pip and the time-translation symmetry of G imply £M{t)^j j d?x'M*{^)S{x'). (19) where g — ~Am{'x)d/c^ , d is the effective nonlinear- ity, equal to half the relevant projection of x''^^ ^^id Afc = kp — ks — ki is the wave-number mismatch. The dimensionless function to(x) describes the distribution of X ■ For example in a periodically-poled material it alternates between ±1. We can take Afc = with- out loss of generality, as the phase oscillation can be in- corporated directly in the envelopes. The propagation directions (±) will be omitted unless needed for clar- ity. Note that for transparent materials x^^-* is real, and X^'^\ujp;ujs+uji) = x^^H^s;^p-Wi) = x^^Hwi;^p-^s)- First-order perturbation theory is sufficient to describe situations in which pairs are produced. For example, if £os , £oi , £op are source- free solutions, then £s = £os+t^s / d'^x'Qs{x;x') xgix')£opix')£lix') + 0{g^). (23) and similar expressions for £i , £p are sufficient to give the lowest-order contribution to the pair-detection rate 14^^^) (X (£l£l£i£sj- Higher-order expansions would be necessary for double-pair production, etc. F. narrow-band frequency Alters Similarly, if a product £i (a;i )£^2 {^2 ) is given by a constant pair source 5^^^ (a:) as fi(a;i)£2(a;2) = / d'^xQi{xi;x)g2{x2;x') x5(2)(x'). (20) then the time-integrated mode-projected component is 1 dti£iMi {tl)£2M2 (^2) d^x'M*{^') PlzP2z J xM;{x')S^^^{x'). (21) E. Coupled wave equations We now introduce a x nonlinearity, which produces a nonlinear polarization that appears as a source term in the propagation equations. We consider three fields, "signal," "idler" and "pump" with carrier frequencies ujs,uJi,ujp and wave- numbers ks,ki,kp, respectively. The respective field envelopes £s,£i,£p evolve according to UpCp Vs£s V,£, ujpg£s£i exp[iAfcz] .2„c- ct.^r^j^fc^i ujfg£p£l exp[— iAfcz] (22) In most down-conversion experiments, some sort of fre- quency filter is used. Assuming this filter is linear and stationary, the field reaching the detector is £'-^Ht) = f dt'F{t - t')£{t') + G{t - t')frcs(t')- (24) Here 5res is a reservoir field required to maintain the field commutation relations. Assuming the reservoir is in the vacuum state, it will not produce detections and can be ignored. Defining HpiU.ts) = (£f'\ti)£i'''\ts)) , the fields that leave the filter obey Hf{U:U) = dt'dt"F,{U-t')Fs{ts-t") X {£,{t')£,{t")) . (25) In the narrowband case, i.e., when the correlation time between signal and idler is much less than the time- scale of the impulse response functions, we can take {£^{t')£s{t")) « A5{t' - t") where the constant A = ]dU{£,{U)£s{ts)). We find HF{U,ts) « A I dt'F,{U-t')Fs{ts-t') = Af{ts-ti). (26) With this, we see that the flux of pairs is W^'\U-U) = ^^,f'"'"° |.4/(t.-t.)P (27) ft2, UJsUJi with a total coincidence rate of A. sum-frequency generation W^(2) ^ dUW(^\t,-U) 2 2 TlsTliC £"q I . 12 ;i2 A\' / dt,|2/(i,-t0|^ hPujs^i \A\^T, eff- (28) We note that Tcff dr2Ts(r2)Ti(-f7) (29) where Ts^i(r2) = | J dt exp[iflt]Fs^i{t)\'^ are the signal and idler filter transmission spectra, respectively. For this reason we refer to Fofr as the effective line-width (in an- gular frequency) for the combined filters. Also important will be the singles rate w^'^ = A;>,,([£f-)(t.)]t£:r^^)(i.)) 2nsCeo dt'dt"F*{ts-t')Fs{ts-t") C I dt'\Fs{t,-t')\^ 2rtsC£o risceo 2fujj, CT, (30) where C = J dt' {£Ut')£s{t")). TeS.s is the effective line- width for the signal filter. IV. RESULTS With the calculational tools described above, we now demonstrate the central results of this paper. We first express the efficiency of continuous-wave sum-frequency generation (SPG) in terms of a mode-overlap integral. This effectively reduces the non-linear optical problem to three uncoupled propagation problems. We then show that the efficiency of parametric down-conversion in the same medium is proportional to the SFG efficiency, for modes with the same shapes but opposite propagation direction. The constant of proportionality is found, al- lowing calculations of absolute efficiency based either on material properties such as x*"^"* or measured SHG effi- ciencies. Similarly, the singles production efficiency is re- lated to difference-frequency generation (DFG) and the collection efficiency is calculated. The same quantities for the degenerate case are also found. We consider first the process of SFG, for un-depleted sig- nal and idler and no input pump. Signal and idler are constant and come from single-modes, £Mp{tp) — — i?j\/;£ Ms C Pz,p X / d^x' M;{x')Tn{x')M^{x')Msix') — —£]\i,£. AluM Mi^Ms c¥.. -ISFG- (31) The conversion efficiency is Pm. Q SFG _ 8L0*npd^\IsFG\' 2^2 2w^d c'eonpUsrii LSFG\ (32) The efficiency of a cw, single-mode source is thus pro- portional to the spatial overlap of the pump, signal, and idler modes, weighted by the nonlinear coupling g. B. non-degenerate parametric down-conversion Next we consider the process of parametric down- conversion. Using Equation (P5)l . we can calculate to first order in g the correlation function {£i{xi)£s{xs)) = ujI j d'^x'Qs{xs,x) X (£oAxi)£l,{x'y^ xg{x')£o^p{x') c^£o J xg,{x,,x')g{x')£o,p{x') (33) For constant pump and single-mode collection we have •^MiMs = j dts {£Mi{ti)£Ms{ts)) ihhJs^id X I d^x'Ml[x')M*{x')fn{x')Mp{x') -£m„Idc- ihiUsLUid (34) We note that Ij^c = ^sfg- Also, the conjugate modes describe backward-propagating fields, as if the source fields were sent through the nonhnear medium in the op- posite direction. Thus if we want to know the brightness of down-conversion when all beams are propagating to the left, it is sufficient to calculate (or measure) the effi- ciency of up-conversion when all beams are propagating to the right. Using equations ((32|) and (|34|) we find 1^ Mi Ah PpQ SFG- (35) generation. We find Pr = c^eonsHinp^ ^ PpPs / d XiS{zi - zq) f3,^, / d^x'g^ix^■,x')mix')Mp{x')M;{x') P P 2luM^ r(«) |2 c-^eonsniUp P'pPsQdfg- (38) C. brightness E. singles rates in PDC We can now consider the brightness of the filtered, single- mode source. The rate of detection of pairs is T^(2) ^ h?U!sUJi ^l^preff — Fpff— — —PpQ H P^SFG (36) This simple expression is the first main result: The rate of pairs is simply the joint collection bandwidth Feff, times the ratio of frequencies, times the pump power, times the up-conversion efficiency QgpQ- Note that the last quan- tity can be calculated if the mode shapes and x''^-'(x) are known, for example in the paper of Boyd and Klein- man, or simulated for more complicated situations. Most importantly, it is directly measurable. D. difference- frequency generation We now consider the classical situation in which pump and signal beam are injected into the crystal and idler is generated. We will see that this directly measurable process is related to the singles generation rate by para- metric down-conversion. The generated idler is E^{x) = ujf j d'^x'Qi{x]x')g(x')Eop{x')£^Xx'). If pump and signal are from modes Mp, Ms, respectively, we find £i{x) = ^^M^{tp)£*MSts) / d'^x'Qi{x]x') xm{x')Mp{x')M*{x'). (37) We can find the rate of detection of singles in the mode Ms by equation ([50]) and using Equation (|A2p C = dts(^£lj^{xs)£Ms{x's] = f dUd^x,d^x',Msixs)M:{x',) x6{z,s - zo)6{z'^ - zo) {£l{xs)£s{x'^)} - I^^^pI''^'- j d^xd^x'Ms[^)g{^)M;{^) \Pz X (fo.(x)4(^')> M;(x').9(x')Mp(x') 2hlOiLUgd i2 / i4 lit/ II \ n . \£mp\ / dV'(5(z"-zo) (3,,, / d3xG,(x";x)M;(x)m(x)Mp(x) (39) so that W^(i) ^ CUJiUJs 32npnsnieo Ppl^DFGi ^ cff,j — r ff p o^''> (40) F. conditional efficiency The conditional efficiency for the idler (probability of col- lecting the idler, given that the signal was collected) is 1^(2 Tcff \Isfg\ W, (1) |7-(s) 12 \-^dfg\ (41) G. degenerate processes The total power generated is Pi = 2cni£o J d^Xi5{zi — Zf))\£i{xi)\'^ where zq indicates a plane downstream of the Up to this point, we have discussed only non-degenerate processes, i.e., those in which the signal and idler fields are distinct and do not interfere. This is always the case for type-II down-conversion, and will be the case for type- I down-conversion if the frequencies and/or directions of propagation are significantly different. We now consider degenerate processes, in which there is only one down- converted field (signal). The above discussion is modified only slightly. The signal and pump evolve by V,£, = io^sg^pSl (42) H. second harmonic generation The calculation of second-harmonic generation (SHG) proceeds exactly as in sum-frequency generation, except for the factor of one half and with all "idler" variables replaced by "signal" variables. Thus we find Pp = P^QSHG where and QsHG ^P^' \r |2 A'-shg\ 2c^eonpn1 IsHG = /'d3xM;(x)m(x)Af,(x)Af,(x). I. average parametric gain (43) (44) (45) term is the phase-independent contribution to the gain, an experimentally accessible quantity. We have = ^^^^^\£0s\'\£p\' J d'xsS{zs - zo) Pz,s / d^x'gs{xs;x') X m(x')Mp(x')M;(x')| c-^n^npEo PqsPpQapg- -\Iapg\ (48) J. degenerate PDC Next we consider the process of degenerate parametric down-conversion, for which £s = £as + ^l / d'^x'Qs{x;x') ^g{x')Eop{x')£l{x'). (49) The other classical process of interest is parametric am- plification of the signal by the pump. The first-order solution for the signal field is £s = Eos+ojI I d'^x'Qs{x;x') xg{x')Sop{x')£l^{x') = ^os+^is. (46) The signal power at the output is r Ps = 2nsceQ / d^XsS{zs ~ ZQ)\£sixs)f = 2nsceo / d^XsS{zs - zq) [\£o,s{xs)\'^ +2Re[£o^s{xs)£ls{xs)] + \£iAxs)\^) ■ (47) The first term is the input signal power Po,s, the second term depends on the relative phase (j)p — 2(j)s, and the last We find the correlation function {£s{xs)£s{x's)) = ojI j d'^x"Qs{xs,x") x(£,A<)£Ux"))g{x") x£o.p(x") (50) at which point it is clear that the only difference from the non-degenerate case of Eq. (|55|) will be the replacement of idler variables with signal variables. We find ^^'^ = y.S-r^PpQsHG = -^PpQsHG- (51) 8 K. Singles rates (degenerate) As before, we can find the rate of detection of singles in the mode Ms by equation (PU)) and C = I dt',(£{,^{Q£Ms{t's)) xS{z'^ - zq)S{z'J - zo) {£l{x',)£s{x'^)) X (fo.(x')4(^")) AC(x").9(x")Af,(x") X /3,,, f d^xGsix";x)M:ix)m{x)Mp{K) The singles rate is thus L. Conditional efficiency (degenerate) (52) (53) where m G {s, i,p}, r is the radial component of x, and q = z — izR where zr is the Rayleigh range, assumed equal for all beams. We assume a periodically-poled ma- terial in which x^^^ (z) alternates with period 2t:/Q so and we approximate m(x) w eyi^[iQz\dcald. From Equation ([5T|1 we find 'SFG ;- TT^ J- L/2 L/2 -iAkz' ^y 2Trr'dr'e^ i^ , ; 2 2i / f^pf^sf^i^j^ J V TT r.L/2 -iAkz' ^y d-' ^ -L/2 iz' — i Zfl) (i?fcZ' + iZR) (56) where Afc = kp — kg — ki — Q and -Rfc = /j_/fc+ and fc± = kp ± (ks + ki). In terms of the dimensionless variables K = AkL, C = z'/L, (^ji = Zfi/L we find 2i IsFG = — a/ttAP^A^TT where T = Ci^ 1/2 dc 'IkC, 27ry_i/2 Mc - *Cfl) (^fcC + <fl) ■ (57) (58) The conditional efhciency is W, (1) -Lcff I^SifGl (54) From Equation (|32p the upconversion efficiency is then n,p Kg Ki Q SFG SttljI c^SonpUsrii kl -Znd^g\T\ (59) A. Boyd and Klienman, 1968 EXAMPLE CALCULATIONS We now illustrate the preceding, general results with a few special cases. We first calculate the overlap integral for co-propagating gaussian beams. This allows us to 1) compare our results to the classical results of Boyd and Kleinman [23|, 2) predict absolute brightness for an im- portant geometry, type-II co-linear down-conversion in quasi-phase-matched material. Also, we compare to a recent calculation of absolute brightness for a specific ge- ometry by Ling, et al. |17| . We consider collinear, frequency-degenerate type-II PDC with circular gaussian beams for signal, idler and pump. We take mode shape functions M,„(x) kmZR 1 TT q ^ik^z^ik„,^^ (55) With this expression we can compare our results to those of Boyd and Kleinman [231 f^'' the case of second- harmonic generation. As that calculation does not in- clude quasi-phase matching, we take Q — 0, and then for any reasonable phase-matching we have Up « rig, kp « 2ks and thus k- = Rk = 0, k+ w 2kp. We note that for i?fc = 0, T becomes equal to the function H of Boyd and Kleinman for zero absorption and walk-off angle. We find Pp = Ank^ c^eon'sUp ZHd'\T\'Pl (60) Boyd and Kleinman find in Equations (2. 16), (2. 17) and (2.20) to (2.24) P2^'-^^d^PiLk,^^\Hf (61) P2 = (62) When converting this expression to MKS units, cP — (P/64TT^eo, and we see that the two calculations agree. B. Type-II collinear brightness Next we make a numerical calculation for frequency- degenerate type-II SPDC, a geometry of current interest for generation of entangled pairs, for example. The integral T must be evaluated numerically. For a 1 cm crystal of PPKTP and a vacuum wavelength A^ = Ai = 800 nm we have {ns,m,np) = (1.844,1.757,1.964) and dcs — 2.4 pm/V so that Rk — 0.04 and the maximum of {zr/L)\T\^ « 0.054 occurs at k « -3.0, Ci?, « 0.18. We find QsFG = 2.0 X 10~^ W~^ Used as a photon-pair source, this same crystal and geometry would yield by Equation VK(2) = FeffFp- :SFG 16 (63) or a pair generation efficiency of Qsfg/^Q = 0.8 pairs (s mW MHz)-^ Note that T^s is the fiher bandwidth in angular frequency. C. Ling, Lamas-Linares, and Kurtsiefer, 2008 Recently, Ling et al. calculated the absolute emission rate into gaussian modes in the thin-crystal limit of (non- periodically-poled) nonlinear material [17| . They arrive to a down-conversion spectral brightness of dR{uJs) _ / rfcffa^a»£^"^(Afc) \ qj^uj (64) du)s \ ^ I 27ms rii where R is the pair collection rate and $(Afc) = I dz f dydxe'^^''Up{r)Us{r)U^{r). (65) Here Um describe the mode shapes of the form Um{^) = f^ik^z^-{x^+y^)/wl ^^^ Q,^ ^ y^2/nWl are normal- ization constants. The field E^ is defined such that i?°P — 2apPp/eQnpC where Pp is the pump power, giving dRJLJs) duj. UJs'^i<P Trc-^eonpUsni Pp \apasa^<^{Ak)\ (66) Assuming the output is collected with narrow-band fil- ters of transmission Tg[u)s),Ti(uJi) for signal and idler, respectively, the integrated rate is R dR{LUs diu. dnTs{ujp/2 + n)Ti{ujp/2-n) (67) where we have assumed dR{ujs)/duJs constant over the width of the filters. For comparison, using Eqs ([32|) and (1361) . we find W-(^) = .^^-^' Pp\Is '^ 2c^eonpnsni FG\ ioff- (68) Then with Equation (|29|) and noting that in the thin crystal limit |/sfg| — \oipCtsOii^{Ak)\, we see that the two results are identical. VI. CONCLUSIONS Using the approach of coupled wave-equations, famil- iar from nonlinear optics, we have calculated the abso- lute brightness and temporal correlations of spontaneous parametric down-conversion in the narrow-band regime. The results are obtained with a Green function method and are generally valid within the paraxial regime. We find that efficiencies of SFG and SPDC can be expressed in terms of mode overlap integrals, and are proportional for corresponding geometries. Also, we find pair time correlations in terms of signal and idler filter impulse re- sponse functions. Results for both degenerate and non- degenerate SPDC are found. Comparison to classical cal- culations by Boyd and Kleinman, and to a recent calcu- lation by Ling et al. show the connection to classical nonlinear optics and "golden rule" -style brightness cal- culations, while considerably generalizing the latter. We expect these results to be important both for designing SPDC sources, as the results of well-known classical cal- culations can be used, and for building and optimizing such sources. APPENDIX A: ALTERNATE PROPAGATOR We can use Equations (|T2l) and p3)l to express the prop- agator as {£{x)£\x')) = \Pt? I d^x"d^x"'5{t" -tf) xSr-tj){£ix")£\x"')) xg*{x"-x)g{x"'-x'). = \A^Pt\^ j d^x"5{t"-tf) y<g*{x"-x)g{x"-x') (Al) 10 Noting that J d-^x"S{t" - tf)g*{x";x)g{x";x') Vgjd'^x"6{z" - zo)g*{x";x)g{x";x') we find {£{x)£Hx')) c^eo d''x"5{z" - zo) xg*{x"-x)g{x";x') (A2) APPENDIX B: LORENTZIAN FILTER A common filter has a Lorentzian transfer function and an exponential impulse response terms, not having a filter there at all), filter becomes fits - U) 2 ti<ts exp[-r,(tj-is)/2] ti >ts (B5) That is, the idler photon will always arrive later, and with a distribution (after the signal arrival) that is precisely the transfer function of the idler-beam filter. Another interesting limit is for matched filters, Fg = F^ = F. Then we find fits - ti) F exp[-F|i«-t,|/2]. (B6) F{t) = -9{t)cM-'Ct/2]. (Bl) The spectral transmission is Tip.) ^ V'^/iV'^ + 4fi2), i.e., unit transmission for constant f , a full- width at half- maximum of AilpwHM = r and an area J dflT{Q) = TrT/2. If we put a filter of this sort in each arm, the output has lits-ti) r^r. dt'eits - t')eiu - 1') fir) X exp[-F,(t, - t')/2] xexphF,(t, -i')/2] TsT^ f exp[-F,T/2] t > 2(F, +F,) lexp[F,T/2] t<0 (B2) The effective bandwidth is r.ff^A fdT\f{T)\ (B3) (B4) It is worth noting that in the limit F^ ^ cxd (the limit of a broad-band filter in the signal beam, or in practical Note that for F^ -^ oo, the detection rate is |^pFi/4, i.e., proportional to the idler filter bandwidth F^. The reverse, s <-> i is also true, of course. From this we can get an idea of the conditional efficiency: The rate for filtered signal with any idler is proportional to F. > F F F., (B7) For example, putting matched filters Fg = F^ = F will give a rate proportional to FiFs/(Fs -I- F^) i.e., half of the rate without the idler filter. This indicates that, of the signal photons that pass the the signal filter, half of their "twin" idler photons do not pass the idler filter. Acknowledgements We thank A. Cere, F. Wolfgramm, G. Molina, A. Haase, and N. Piro for helpful discussions. This work was sup- ported by the Spanish MEC under the ILUMA project (Ref. FIS2008-01051), the Consohder-Ingenio 2010 Project "QOIT" and by Marie Curie RTN "EMALI." [1] M. Fejer, G. Magel, D. Jundt, and R. Byer, IEEE JOUR- [7] NAL OF QUANTUM ELECTRONICS 28, 2631 (1992). [2] P. Kwiat, K. Mattle, H. Weinfurter, A. Zeilinger, [8] A. Sorgicnko, and Y. Shili, PHYSICAL REVIEW LET- TERS 75, 4337 (1995). [9] [3] P. Kwiat, E. Waks, A. White, I. Appelbaum, and P. Eberhard, PHYSICAL REVIEW A 60, R773 (1999). [4] G. Giorgi, G. Di nepi, P. Mataloni, and F. De martini, [10] LASER PHYSICS 13, 350 (2003). [5] M. Fiorentino, G. Messin, C. Kuklcwicz, F. Wong, and [11] J. Sliapiro, PHYSICAL REVIEW A 69, 041801 (2004). [6] M. Pelton, P. Marsden, D. Ljunggren, M. Tengner, A. Karlsson, A. Fragemann, C. Canalias, and F. Laurell, [12] OPTICS EXPRESS 12, 3573 (2004). C. Kuklewicz, M. Fiorentino, G. Messin, F. Wong, and J. Sliapiro, PHYSICAL REVIEW A 69, 013807 (2004). M. Fiorentino, C. Kuklewicz, and F. Wong, OPTICS EX- PRESS 13, 127 (2005). F. Wolfgramm, X. Xing, A. Cere, A. Predojevic, A. Stein- berg, and M. Mitchell, OPTICS EXPRESS 16, 18145 (2008). J. O'brien, G. Pryde, A. White, T. Ralph, and D. Bran- ning. Nature 426, 264 (2003). S. Groblacher, T. Paterek, R. Kaltenbaek, C. Brukner, M. Zukowski, M. Aspelmeyer, and A. Zeilinger, Nature 446, 871 (2007). B. Higgins, D. Berry, S. Bartlett, H. Wiseman, and G. Pryde, Nature 450, 393 (2007). 11 [13] S. Sauge, M. Swillo, S. Albert-Seifried, G. Xavier, [23] J. Waldeback, M. Tengner, D. Ljunggren, and A. Karls- [24] son, OPTICS EXPRESS 15, 6926 (2007). [14] A. Haase, N. Piro, J. Eschner, and M. W. Mitchell, OP- [25] TICS LETTERS 34, 5557 (2009). [15] D. Ljunggren and M. Tengner, PHYSICAL REVIEW A [26] 72, 062301 (2005). [16] C. Kurtsiefer, M. Oberparleiter, and H. Weinfurter, [27] PHYSICAL REVIEW A 64, art. no. (2001). [17] A. Ling, A. Lamas-Linares, and C. Kurtsiefer, PHYSI- CAL REVIEW A 77, 043834 (2008). [28] [18] D. Kleinman, PHYSICAL REVIEW 174, 1027 (1968). [19] Y. Zeldovich and D. Klyshko, JETP LETTERS-USSR 9, [29] 40 (1969). [20] Y. R. Shen, The Principles of nonlinear optics (Wiley- [30] Interscience, New York, 1984), 1st ed. [21] D. Klyshko, ZHURNAL EKSPERIMENTALNOI I TEO- RETICHESKOI FIZIKI 94, 82 (1988). [22] D. Burnham and D. Weinberg, PHYSICAL REVIEW LETTERS 25, 84 (1970). B. MoUow, PHYSICAL REVIEW A 8, 2684 (1973). C. Hong and L. Mandel, PHYSICAL REVIEW A 31, 2409 (1985). R. Ghosh, C. Hong, Z. Ou, and L. Mandel, PHYSICAL REVIEW A 34, 3962 (1986). A. Valencia, A. Cere, X. Shi, G. Molina- Terriza, and J. Torres, PHYSICAL REVIEW LETTERS 99 (2007). P. Mosley, J. Lundeen, B. Smith, P. Wasylczyk, A. U'ren, C. Silberhorn, and I. Walmsley, PHYSICAL REVIEW LETTERS 100, 133601 (2008). R. W. Boyd, Nonlinear Optics (Academic Press, Burling- ton, MA, 2008), 3rd ed. G. Boyd and D. Kleinman, JOURNAL OF APPLIED PHYSICS 39, 3597 (1968). We note that Hamiltonian-based treatments often do not agree about the constant preceding the integral J d^xx E^, not even its sign. This question does not present a problem for the present calculation, which does not employ a Hamiltonian.